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The study of internuclear potentials for heavy-ion 
collisions is of fundamental importance for the forma- 
tion of superheavy elements and nuclei far from stabil- 
ity. While asymptotically such potentials are determined 
from Coulomb and centrifugal interactions, the short dis- 
tance behavior strongly depends on the nuclear surface 
properties and the readjustments of the combined nuclear 
system, resulting in potential pockets, which determine 
the characteristics of the compound nuclear system. 

Among the various approaches for calculating ion-ion 
potentials are: 1} Phenomenological models such as the 
Bass model pj 0, the proximity potential 0, S IE BLj 
and potentials obtained via the double-folding method 0, 
IE EH ED ■ Some of these potentials have been fitted to 
experimental fusion barrier heights and have been re- 
markably successful in describing scattering data. 2) 
Semi-microscopic and full microscopic calculations such 
as the macroscopic-microscopic method 0, 0, fl3| . 
the asymmetric two-center shell-model |14|. constrained 
Hartree-Fock (CHF) with a constraint on the quadrupole 
moment or some other definition of the internuclear 
distance [la. Ha . and other mean-field based calcula- 
tions [17 I IT a ligf . 

One common physical assumption used in many of the 
semi-microscopic calculations is the use of the frozen den- 
sity or the sudden approximation. As the name suggests, 
in this approximation the nuclear densities are unchanged 
during the computation of the ion-ion potential as a func- 
tion of the internuclear distance. On the other hand, the 
microscopic calculations follow a minimum energy path 
and allow for the rearrangement of the nuclear densities 
as the relevant collective parameter changes. As it was 
pointed out in Ref. CHF calculations seldom pro- 
duce the correct saddle-point since the system can follow 
any one of the minimum potential valleys in the multi- 
dimensional potential energy surface. In this paper, we 
shall call this the static adiabatic approximation since a 
real adiabatic calculation would involve a fully dynamical 
calculation, thus also including the effects of dynamical 
rearrangements. 

One conclusion that may be reached from the discus- 
sion above is that ultimately we would like to have an 
approach for calculating internuclear potentials which is 
time-dependent and is unrestricted in the choice of collec- 
tive variables. In this paper we provide such an approach 
in which time-dependent Hartree-Fock (TDHF) is used 



for the nuclear dynamics and the potential energy is cal- 
culated by constraining the time-dependent density. 

The density constraint is a novel numerical method 
that was developed in the mid 1980's [2CJ, |2l| and was 
used to provide a microscopic description of the for- 
mation of shape resonances in light systems [2l|. In 
this approach the TDHF time-evolution takes place with 
no restrictions. At certain times during the evolution 
the instantaneous density is used to perform a static 
Hartree-Fock minimization while holding the total den- 
sity constrained to be the instantaneous TDHF density. 
In essence, this provides us with the TDHF dynamical 
path in relation to the multi-dimensional static energy 
surface of the combined nuclear system. Since we are 
constraining the total density all moments are simulta- 
neously constrained. In the traditional CHF notation 
this corresponds to the replacement 

XQ — > A/5 . (1) 

The numerical procedure for implementing this con- 
straint and the method for steerin g t he solution to 
Ptdhf (r, t) is discussed in Refs. pfl l2lj. The conver- 
gence property is as good if not better than in the tra- 
ditional CHF calculations with a constraint on a single 
collective degree of freedom. 

In this paper, we shall call the energy of the system ob- 
tained by the density constraint method, Edc(R), where 
the dependence is on the instantaneous internuclear sepa- 
ration, R(t). Since this quantity contains no translational 
kinetic energy (taken out by the static minimization) it 
is actually a potential as we show below. We define the 
excitation energy of the system as 

E (R) = -Btdhf — Tr — Et)c(R) , (2) 
where £tdhf is the total TDHF energy 

£tdhf = f d 3 r H(r,t) , (3) 

which is conserved throughout the calculation, and Tr is 
the instantaneous translational energy between the two 
nuclei 

Tr = \^R 2 , (4) 

with \i being the reduced mass of the system, and R is 
the velocity associated with the internuclear separation 
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coordinate R(t). At the same time the total TDHF en- 
ergy can be written in terms of the excitation energy as 



E TunF =T R + V(R)+E*(R) , 



(5) 



such that when combined with Eq. @ shows that 
V(R) = Evc(R)- However, the density constrained po- 
tential still contains the binding energies of individual 
nuclei, which should be subtracted out; 



V(R) - V(R) = E DC (R) - E Al - E A2 



(6) 



Eq. © is the internuclear potential and contains no free 
parameters. Given an effective nuclear interaction, such 
as the Skyrme force, V(R) can be constructed by per- 
forming a TDHF evolution and minimizing the energy 
at certain times to obtain Edc(R) 7 while E Al and E^ 2 
are the results of a static Hartree-Fock calculation with 
the same effective interaction. One can see that the ex- 
pression also has the correct asymptotic behavior since 
numerically for large R we exactly get 



max ! 



E Al + E A2 + 



Z±Z 2 e 2 

^max 



such that 



V(R ma , x ) 



Z!Z 2 e 2 



(7) 



(8) 



so that normalization of V(R) is not necessary. 
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Also shown in Fig. are two widely used phenomeno- 
logical potentials, the standard proximity potential for 
two spherical nuclei 0, 0, Hi an d the double- folding po- 
tential with M3Y effective NN interaction 0, H |gLflo| . 
We evaluate the double-folding integral for the stron g nu - 
clear and Coulomb interaction in momentum space |10| . 
For the charge and matter densities we utilized gener- 
alized Fermi distributions whose parameters were deter- 
mined from electron scattering experiments [25|. The 
double-folding potential agrees almost perfectly with the 
DC-TDHF approach for distances R > 6 fm. At smaller 
distances, the double-folding potential tends to overesti- 
mate the nuclear interaction as a result of the (unphys- 
ical) frozen density approximation. The classical TDHF 
curve corresponds to the older definition 



V(R) = E c 



(9) 



used in Ref. [24j. 

We have also repeated the above calculations for differ- 
ent center-of-mass energies. In this case, we find that the 
results at the barrier do not appreciably change but the 
depth of the potential increases for lower energies. For 
Ec.m. — 12 MeV the potential is about 1.5 MeV deeper 
than the one at 34 MeV. One comment is required re- 
garding the calculation of the internuclear separation R. 
As usual, this quantity becomes somewhat unclear for a 
strongly overlapping system. In our case we use the stan- 
dard TDHF approach of finding left and right dividing 
planes and computing the centers of the density in these 
two halves and thus the separation. 
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FIG. 1: Internuclear potential obtained from Eq. |S| and 

various model calculations for the head-on collision of the 
i6q + i6q system at Ec m _ 34 MeV 



FIG. 2: Internuclear potential obtained from Eq.^shown for 
the entire time evolution of the 16 0+ 16 system at -E c .m. = 
34 MeV. 



We have carried out a number of TDHF calculations 
with accompanying density constraint calculations to 
compute V(R) given by Eq. |JB}. A detailed description 
of our new three-dimensional unrestricted TDHF code 
has recently been published in Ref [2^]. For the effec- 
tive interaction we have used the Skyrme SLy5 force j2^| 
including all of the time-odd terms. In Fig. ^ we show 
the result of our calculation for the head-on (zero impact 
parameter) collision of 16 0+ 16 at E c . m . = 34 MeV. 



In TDHF, fusion occurs when the relative kinetic en- 
ergy in the entrance channel is entirely converted into 
internal excitations of a single well defined compound 
nucleus. The dissipation of the relative kinetic energy 
into internal excitations is due to the collisions of the nu- 
cleons with the "walls" of the self-consistent mean-field 
potential. TDHF studies demonstrate that the random- 
ization of the single-particle motion occurs through re- 
peated exchange of nucleons from one nucleus into the 
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other. Consequently, the equilibration of excitations is 
very slow and it is sensitive to the details of the evolution 
of the shape of the composite system. This is in contrast 
to most classical pictures of nuclear fusion, which gener- 
ally assume near instantaneous, isotropic equilibration. 
This equilibration can be observed in the DC-TDHF ap- 
proach by tracking the evolution of the excitation energy 
in time for a system on the way to fusion, or alternately 
one can examine the change in the internuclear potential 
for the compound system. In Fig. |3 we show this for 
the 16 0+ 16 system corresponding to the case shown 
in Fig. ^ After passing the first minimum the system 
falls back to a second minimum and climbs up the po- 
tential barrier, but it cannot overcome the barrier due 
to some of the energy being converted into internal ex- 
citations, consequently it falls back to a third minimum, 
fourth minimum, and so on, until complete equilibration. 
This potential ladder effect is characteristic for all fusing 
systems. 
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FIG. 3: Internuclear potential obtained from Eq.|S|shown for 
the evolution of the 16 0+ 22 Ne system at E c . m . = 50 MeV. 
The two curves correspond to different orientations of the Ne 
nucleus. 



Finally, we have performed calculations for the 



16 0+ 22 Ne system at E c . m . = 50 MeV. In this case the 
22 Ne nucleus shows a strong axial deformation, which 
can have different orientations with respect to the col- 
lision axis. We have recently reported a procedure for 
doing such calculations within the TDHF framework in 
Ref. |2(j ■ Here, we show the result of our potential calcu- 
lations for two orientations of the 22 Ne nucleus, one for 
which the symmetry axis of the Ne is aligned with the 
collision axis and the other for which the symmetry axis 
is perpendicular to the collision axis. In TDHF, for the 
first case (aligned with collision axis) we see no fusion at 
this energy, whereas for the perpendicular alignment the 
system fuses. In Fig.[3]we show the results of our poten- 
tial calculations. We observe that the barrier height and 
its position in R space, and the potential minimum and 
its position are considerably different for the two orien- 
tations. 

In summary, we have presented a method for calculat- 
ing internuclear potentials directly from the TDHF time- 
evolution of the colliding system. The method uses the 
density constraint to trace the TDHF trajectory in rela- 
tion to the static multi-dimensional energy surface of the 
combined system. Since the TDHF evolution is unhin- 
dered all of the collective dynamics associated with the 
evolution are included in the calculations. We believe 
this provides a unique way to calculate ion-ion potentials 
from mean-field calculations. Of course, we can only per- 
form such calculations for energies above the barrier since 
TDHF is a semi-classical theory in this regard. Finally, 
the results are expected to be only as good as the TDHF 
description of the particular system under study. 

This work has been supported by the U.S. Depart- 
ment of Energy under grant No. DE-FG02-96ER40963 
with Vanderbilt University. Some of the numerical cal- 
culations were carried out at the IBM-RS/6000 SP su- 
percomputer of the National Energy Research Scientific 
Computing Center which is supported by the Office of 
Science of the U.S. Department of Energy. 



[1] R. Bass, Nucl. Phys. A231, 45 (1974). 

[2] R. Bass, Nuclear Reactions with Heavy Ions, (Springer- 

Verlag, New York, 1980). 
[3] J. Blocki, J. Randrup, W. J. Swiatecki, and C. F. Tsang, 

Ann. Phys. (N.Y.) 105, 427 (1977). 
[4] J. Randrup and J. S. Vaagen, Phys. Lett. B77, 170 

(1978) . 

[5] M. Seiwert, W. Greiner, V. E. Oberacker. and M. J. 

Rhoades-Brown, Phys. Rev. C 29, 477 (1984). 
[6] J. R. Birkelund and J. R. Huizenga, Phys. Rev. C 17, 

126 (1978). 

[7] G. R. Satchler and W. G. Love, Phys. Rep. 55, 183 

(1979) . 

[8] G. Bertsch, J. Borysowicz, H. McManus, and W. G. Love, 
Nucl. Phys. A284, 399 (1977). 



[9] M. J. Rhoades-Brown and V. E. Oberacker, Phys. Rev. 

Lett. 50, 1435 (1983). 
[10] M. J. Rhoades-Brown, V. E. Oberacker, M. Seiwert, and 

W. Greiner, Z. Phys. A 310, 287 (1983). 
[11] H. J. Krappe, J. R. Nix, and A. J. Sierk, Phys. Rev. C 

20, 992 (1979). 

[12] P. Moller, A. J. Sierk, and A. Iwamoto, Phys. Rev. Lett. 
92, 072501 (2004). 

[13] T. Ichikawa, A. Iwamoto, P. Moller, and A. J. Sierk, 
Phys. Rev. C 71, 044608 (2005). 

[14] J. Maruhn and W. Greiner, Z. Phys. 251, 431 (1972). 

[15] P. G. Zint and U. Mosel, Phys. Rev. C 14, 1488 (1976). 

[16] T. Biirvenich, M. Bender, J. A. Maruhn, and P.-G. Rein- 
hard, Phys. Rev. C 69, 014307 (2004). 

[17] V. Yu. Denisov and W. Norenberg, Eur. Phys. J. A 15, 



4 



375 (2002). 

[18] A. Dobrowolski, K. Pomorski, and J. Bartel, Nuclear 
Physics A729, 713 (2003). 

[19] H. Goutte, J. F. Berger, P. Casoli, and D. Gogny, Phys. 
Rev. C 71, 024316 (2005). 

[20] R. Y. Cusson, P. -G. Reinhard, M. R. Strayer, J. A. 
Maruhn, and W. Greiner, Z. Phys. A 320, 475 (1985). 

[21] A. S. Umar, M. R. Strayer, R. Y. Cusson, P. -G. Rein- 
hard, and D. A. Bromley, Phys. Rev. C 32, 172 (1985). 

[22] A. S. Umar and V. E. Oberacker, Phys. Rev. C 73, 
054607 (2006). 



[23] E. Chabanat, P. Bonche, P. Haensel, J. Meyer and R. 

Schaeffer, Nucl. Phys. A635, 231 (1998); Nucl. Phys. 

A643, 441 (1998). 
[24] S. E. Koonin, K. T. R. Davies, V. Maruhn-Rezwani, H. 

Feldmeier, S. J. Krieger, and J. W. Negele, Phys. Rev. C 

15, 1359 (1977). 
[25] I. Sick and J. S. McCarthy, Nucl. Phys. A150, 631 

(1970). 

[26] A. S. Umar and V. E. Oberacker, (submitted to Phys. 
Rev. C) (2006);pittp://arxiv.org/abs/nucl-th/0604010 . 



